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(57) Abstract: The invention relates to a method for geothechnically characterising an underwater floor (B) involving at least one 
phase, wherein an acoustic wave is emitted by a source (1) and detected by a receiver streamer (2) and an analysis phase, wherein 
signals produced by receivers (21 -24) are used for determining the physical parameters of a strata forming the floor (B). The inventive 
method uses the angle and frequency-dependence of the acoustic wave reflection coefficient on the interface formed by each stratum. 



(57) Abrege : L invention conceme un procede de caracterisation geotechnique d'un fond sub-aquatique (B) comprenant au moins 
une phase d' investigation dans laquelle une onde acoustique est emise par une source (1) et detectee par une flute (2) de recepteurs, 
et une phase d' analyse dans laquelle les signaux produits par les recepteurs (21a 24) sont exploites pour determiner des parametres 
physiques des strates constituant le fond (B). Le procede de Tinvention exploite la dependance angulaire et frequentielle du coeffi- 
cient de reflexion de Tonde acoustique sur T interface que forme chacune des strates. 
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PROCEDE DE CAH2KCTERISATION 6EOTECHNIQUE D'UN FOND 
SUBAQUATIQUE , UTIIiISANT UNE ONDE ACOUSTIQUE A LARGE BANDE 

ET A RASANCE MULTIPLE. 

5 L' invention concerne, de fagon generale^. les techniques 
de prospection par excitation sonique. 

Plus precisement, 1' invention concerne un procede de 
caracterisation geotechnique d'un fond sub-aquatique, tel 

10 qu'un fond sous-marin^ recouvert d'une nappe d'eau de 
hauteur tot ale determinee, et comprenant une pluralite de 
strates formant entre elles des interfaces de separation, 
presentant chacune des parametres physiques propres et 
s ' etendant a des profondeurs diverses en dessous d'une 

15 premiere strate qui forme une interface de reparation 
avec I'eau, ce procede comprenant au moins une phase 
d' investigation, une phase de pre-traitement et une phase 
d " analyse, la phase d' invest iga:t ion incluant elle-meme 
une operation d' excitation acoustique, mise en oeuvre en 

20 emettant au moins vers le fond, depuis une source 
acoustique immergee, une onde acoustique incidente de 
signature f requentielle connue, une operation 

d * acquisition mise en oeuvre en produisant, au moyen d'une 
flute d'au moins quatre recepteurs immerges, des signaux 

25 de mesure respectifs resultant d'une detection d'ondes 
acoustiques respectives reflechies par le fond, et une 
operation de releve de donnees geometriques, la flute 
etant au moins approximativement alignee sur la source et 
distante du fond, les recepteurs etant espaces les uns 

30 des autres, et I'onde incidente etant emise dans chaque 
fenetre d'une succession de fen§tres temporelles 
d' emission disjointes, et presentant une variation 
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temporelle de frequence a 1 " interieur de chaque fehStre 
d' emission. 

Un precede de ce type est par exemple decrit dans le 
5 brevet US 6 151 556. 

D'autres precedes sent egalement connus de 1 ' homme du 
metier, et notamment ceux que decrivent les documents de 
brevets EP - 0 215 703 et EP - 0 553 053. 

10 

De fagon plus g6nerale, les systemes acoustiques connus 
pour 1' investigation des fonds sont soit des systemes 
d^imagerie, soit des systemes sismiques. 

15 Les systemes d^imagerie classiques utilisent des ondes 
acoustiques de frequences elevees, . typiquement 
super ieures a 30kHz, emises vers le fond et permettant 
d^etablir une image de la surface du fond ind6pendamment 
du sous-sol en presence en se basant sur le mecanisme de 

20 retro-diffusion . 

v.* 

Les systemes sismiques usuels du type precedemment defini 
utilisent, en tant qu'emetteur, un boomer produisant un 
signal acoustique a tres basse frequence et presentant 
25 done I'avantage de penetrer profondement dans le sous- 
sol. 

Le signal issu des multiples reflexions et recueilli sur 
une flute d' hydrophones et/ou de geophones tractee par un 
30 navire en synchronisme avec le boomer est ainsi 
representatif des couches prof ondes du sous-sol. 
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En revanche^ les caracteristiques f requentielles du 
signal emis a tres basse frequence sont mal contr61ees, 
de sorte que ces systemes ne fournissent que des 
indications qualitatives des strates profondes du sous- 
5 sol^ a des profondeurs typiquement superieures a 10 
metres, et avec des resolutions verticales typiquement 
superieures a 1 metre. 

Les sondeurs de sediment recents presentent I'avantage 
10 d'utiliser des signaux controles ("chirp" ou sources 
parametriques) qui insonifient le fond verticalement a 
forte puissance. Toutefois, la reception des signaux, 
assures par le mSme transducteur , limite 
signif icativement les performances, pour des raisons 
15 instrumentales (saturation electronique) et physiques 
(melange de phenomenes de retro-diffusion et de reflexion 
speculaire, ambigulte epaisseur/vitesse dans les 
couches) . 

20 Dans ce contexte, la presente invention vise a proposer 
un precede repondant au besoin, aujourd'hui non 
satisfait,- de permettre la determination quantifiee des 
proprietes geotechniques du sous-sol d'un fond sub- 
aquatique sur les premiers metres, typiquement sur les 

25 cinq premiers metres au moins, et avec une resolution au 
plus egale a cinquante centimetres. 

A cette fin, le precede de 1' invention, par ailleurs 
conforme a la definition generique qu'en donne le 
30 preambule ci-dessus, est essentiellement caracterise en 
ce que la phase d" investigation inclut une operation de 
deplacement, concomitante avec les operations 
d' excitation et d' acquisition, et mise en oeuvre en 
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deplagant simultanement la source acoustique et la flute 
dans la nappe d'eau, en ce que la phase de pre-traitement 
est mise en oeuvre en deduisant, a partir des signaux de 
mesure, de la signature de 1 ' onde acoustique incidente, 
5 et des donnees geometriques, des signaux de reflexion 
ayant atteint des recepteurs differents a la suite d'une 
meme emission d'onde incidente, et en regroupant, sous 
forme de series de signaux de reflexion recales, des 
signaux de reflexion ayant atteint, sous des angles de 

10 rasance differents, des recepteurs differents en 
provenance d'une mSme zone du fond a la suite de 
differentes emissions de I'onde incidente, et en ce que 
la phase d' analyse est mise en oeuvre en deduisant les 
parametres physiques de chaque strate du fond, a partir 

15 des series de signaux de reflexion r^cales, par 
exploitation de la dependance angulaire et f requentielle 
du coefficient de reflexion de I'onde incidente sur 
chaque interface, chaque signal de reflexion recale etant 
interprets comme une somme d'arrivees retardees et 

20 attenuees de rayons propres correspondant a la reflexion 
de I'onde incidente par les differentes strates de la 
zone du fond a laquelle correspond ce signal de reflexion 
recale. 

25 En pratique, chaque onde reflechie atteignant le centre 
de la flute dans une fenetre temporelle. de reception est 
avantageusement issue de la reflexion d'une onde 
incidente ayant atteint le fond sous un angle de rasance 
au plus egal a 70 degres et au moins egal A 10 degres, et 

30 de preference compris entre 25 et 65 degres . 

Chaque recepteur est par exemple separe de . la source 
acoustique par une distance comprise entre environ 1 fois 
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et environ 4 fois la hauteur d'' eau moyenne sous les 
instruments^ chaque recepteur pouvant par ailleurs etre 
espace d'un recepteur voisin par une distance comprise 
entre environ 0.375 fois et environ 1 fois la hauteur 
5 moyenne suivant le nombre de ces recepteurs. 

Selon un autre mode d' evaluation des distances, chaque 
recepteur est typiquement separe d'un recepteur voisin 
par une distance correspondant a une difference d' angle 
10 de rasance au moins egale a 2 degres, et de preference au 
moins egale a 5 degres, des signaux de reflexion 
respectivement regus par ce recepteur et par le recepteur 
voisin. 

15 De preference, 1 ' onde incidente presente, dans, chaque 
fen^tre d' emission, au moins une composante de frequence 
relativement basse, dont la frequence est au moins egale 
a 100 Hz, cette onde incidente presentant utilement 
aussi, dans chaque fenetre d' emission, au moins une 

20 composante de frequence relativement haute^ dont la 
frequence est au plus egale a 8 kHz. 

. Dans un mode de realisation possible de 1' invention, 
I'onde incidente peut ainsi presenter, • dans chaque 
25 fenetre d' emission, une frequence continument variable 
entre une composante de frequence relativement basse et 
une composante de frequence relativement haute. 

Pour bbtenir a la fois une information suffisamment riche 
30 et une separation satisf aisante des ondes reflechies, les 
fenetres temporelles d^emission ont chacune par -'exemple 
une duree au moins egale a 0.1 seconde, et de preference 
au moins egale a 0.5 seconde, et sont repetees 
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periodiquement avec une periode par example au moins 
egale a 0.5 seconde et de preference au moins egale a 1 
seconde, 

5 Dans ' son mode de realisation prefere, le procede de 
1 ' invention comprend, en amont de la phase d' analyse^ une 
operation de determination de la profondeur d' immersion 
de la source acoustique et/ou de la profondeur 
d' immersion d^un point de la flute et/ou de la hauteur 
10 d'eau. 

La phase d' analyse est par exemple mise en cBuvre en 
parcourant, de fagon iterative^ une boucle incluant une 
operation de mod61isation au cours de laquelle sont 

15 calcules, a partir de parametres physiques fictifs 
virtuellement attribues a la premiere strate au moins de 
la • pluralite de strates, et pour les differents 
recepteurs^. des signaux de reflexion recales fictifs 
respectifs issus de signaux de mesure virtuellement 

20 produits par ces recepteurs en presence de ces^. parametres 
fictifs une operation de soustraction au cours de 
laquelle sont formes des ecarts constates entre les 
signaux de reflexion recales respectifs issus des 
signaux de mesure reellement produits par ces recepteurs 

25 et les signaux de reflexion recales fictifs 
correspondants, et une operation d' optimisation et de 
retroaction a 1' issue de laquelle les parametres 
physiques fictifs sont' revises en fonction des ecarts 
constates . 

30 

Le procede de 1' invention permet ains-l une 
caracterisation geotechnique de chaque strate d'une 
pluralite de strates d'un fond sub-aquatique au moyen de 
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parametres physiques choisis dans 1' ensemble de 
paraitietres qui comprend I'epaisseur de cette strate, la 
densite de cette strate, la vitesse du son dans cette 
strate, 1 ' attenuation du son dans cette strate, la 
5 rugosite de cette strate^ et un gradient de celerite du 
son a I'interieur de cette strate. 

Dans son mode de realisation prefere, le precede de 
1' invention permet d' identifier les differentes strates 
10 du sous-sol sub-aquatique avec des resolutions verticales 
de I'ordre de dix centimetres . 

D'autres caracteristiques ' et avantages de 1' invention 
ressortiront clairement, de la description qui en est 
15 faite ci-apres, a titre indicatif et nullement limitatif, 
en reference aux dessins annexes^ dans lesquels : 

- la figure 1 est une vue schematique en perspective d'un 
equipement mettant en oeuvre le proced6 de 1 ' invention ; - 

20 

- la figure 2 est un diagrairane schematique presentant de 
fagon approximativement chronologique les principales 
phases et operations du precede de 1' invention ; 

25 - la figure 3^ composee des' schemas 3A a 3C, illustre de 
maniere intuitive et partielle I'etape REGAL de recalage 
et de regroupement de signaux de reflexion mise en oeuvre 
pendant le deroulement du precede de 1' invention, cette 
etape conduisant notarament a regrouper des signaux de 

30 reflexion provenant d'une m§me zone du fond; 

- la figure 4, composee des schemas 4A et 4B, complete 
1' illustration de I'' etape REGAL; 
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- la figure 5;. cortiposee des schemas 5A et 5B, illustre 
partiellement les mecanismes de propagation acoustique 
multiple dans I'eau, et leur exploitation dans 

5 1' invention;. 

- la figure 6 illustre partiellement les mecanismes de 
propagation acoustique multiple dans le fond sub- 
aquatique; 

10 

- les figures 7A a 7 J sont des diagrammes representant , 
en fonction du temps t exprime en nombre d'' echantillons, 
1' amplitude relative A de signaux de mesure 
respectivement produits par 10 recepteurs successifs de 

15 plus en plus eloignes de la source acoustique et separes 
les uns des autres par une distance importante^ et 
mettant en evidence 1' influence des angles avec lesquels 
les ondes se ref lechissent au fond ; et 

20 - les figures 8A a 8 J sont des diagrammes representant 
en fonction du temps t exprime en nombre d'' echantillons^. 
1' amplitude relative A- de signaux de mesure 
respectivement produits par 5 recepteurs successifs de 
plus en plus eloignes de la source acoustique et pour une 

25 modulation de basses frequences a gauche et pour une 
modulation de frequences hautes a droite^ mettant en 
evidence 1'' influence des frequences utilisees. 

Notations 

30 Le lecteur trouvera ci-apres quelques-uns des symboles 
utilises pour identifier les differentes ''entites 
impliquees dans la mise en oeuvre de 1' invention : 
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"i" est un indice compteur utilise pour numeroter chaque 
hydrophone en le distinguant ainsi des autres 
hydrophones; 

5 "j" est un indice compteur utilise pour numeroter chaque 
strate du fond sub-aquatique en la distinguant ainsi des 
autres strates du fond; 

"J" est un nombre entier designant le nombre total de 
10 strates observables; 

"k" est un indice compteur utilise pour numeroter chaque 
rayon propre d'une onde acoustique en le distinguant des 
autres rayons propres; 

•15 

"1" est un indice compteur utilise pour numeroter chaque 
emission d'une onde acoustique incidente (encore appelee 
"tir") en la distinguant ainsi des autres emissions d'une 
onde acoustique; 

20 

"n" est un indice compteur utilise pour numeroter chaque 
boucle d' iteration de la phase d' analyse en la 
distinguant ainsi des autres boucles d' iteration; 

25 "e" est un indice attribue a une valeur finale exacte 
d'lin parametre; 

s (t) designe 1 ' onde acoustique incidente; 



Si(t,l) designe le signal de mesure regu sur le i 
30 hydrophone a la suite du i^^^^ tir; 
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ei(t/l) designs I'enveloppe de 1 ' intercorrelation entre 
I'onde acoustique incidente s(t) et le signal de mesure 
Si(t,l); 

5 rfi(t,l) designe le signal de reflexion associe au signal 
de mesure Si(t,l) par application a I'enveloppe ei(t,l) 
d'une fenetre de reception; 

Rfi(t,B) designe un signal de reflexion recale;. issu du 
10 signal de reflexion rfi(t^l) par une .etape de recalage 
REGAL et associe au fond B; 

RfSi(t, B) designe un signal de reflexion recale fictif^ 
obtenu par simulation; 

15 

Rfsi^(t^ B) designe le k^^"^® rayon propre du signal de 
reflexion recale Rfi(t,B); et 

designe 1' angle de rasance sur la j^^"^® strate du 
20 k^^K^e rayon propre du signal de reflexion recuei,illi sur le 
j^^me hydrophone. 

Comme indique precedeminent , 1' invention concerne (figure 
1) un precede de caracterisation geotechnique d'un fond 
25 sub-aquatique B^ tel qu'un fond sous-marin;. recouvert 
d'une nappe d'eau de hauteur totale H determinee. 

Ce fond B comprend une pluralite de strates telles que 
Fl;. F2y etc. (figure 6) ^ generiquement notees Fj et 
30 presentant chacune des parametres physiques propre^. 
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La premiere strate Fl forme une interface avec I'eau, les 
autres strates s ' etendant en dessous de cette premiere 
strate a des profondeurs diverses. 

5 Le precede de 1' invention vise, pour chacune des strates 
Fj d'un ensemble de strates comprenant typiquement, mais 
non limitativement , entre deux et cinq strates, a 
determiner plusieurs parametres distinctifs de cette 
strate, tels que I'epaisseur hj de cette strate, sa 

10 densite pj, la vitesse du son Cj dans cette strate, et 
plus pr^cisement la vitesse Cpj pour des ondes de 
compression et la vitesse Csj pour des ondes de 
cisaillement , 1 ' attenuation du son aj dans cette strate, 
et plus precisement 1 ' attenuation apj pour des ondes de 

15 compression et 1 ' attenuation asj pour des ondes de 
cisaillement, ainsi, eventuellement , que la rugosit6 dans 
cette strate, ou un gradient de celerite du son. 

Le precede de 1' invention comprend essentiellement une 
20 phase d' investigation INVESTIG, une phase de pre- 
traitements PRE_TRAIT et une phase d' analyse ANA. 

La phase d' investigation INVESTIG inclut elle-meme une 
operation de deplacement DEPL, une operation d' excitation 
25 acoustique EXCIT, une operation d' acquisition ACQUI, et 
une operation RELEV de releve de donnees non acoustiques, 
ces operations etant mises en oeuvre de fagon globalement 
simultanee - 

L^ operation de deplacement DEPL consiste a deplacer dans 
30 la nappe d'eau, au moyen d'un navire N et a vitesse 
r6duite idealement constante, une source acoustique 1 et 
une flute 2 au-dessus du fond B, la flute 2 comprenant 
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par exemple trois, quatre, ou cinq recepteurs acoustiques 
21 a 24, encore appeles hydrophones . 

L' operation RELEV de releve de donnees non acoustiques 
5 comprend elle-ineme une operation CAPT de capture de 
donnees environnementales et une operation de reperage 
REFER. 

L'' operation CAPT de capture de donnees environnementales 
10 consiste a enregistrer a chaque instant t les profondeurs 
de la source acoustique 1 et d' un point de la flute 2, 
ainsi que la hauteur d' eau H. 

operation de reperage REFER consiste a enregistrer, a 
15 chaque instant t, les coordonnees X et Y du point du 
navire N. 

operation EXCIT d' excitation acoustique est mise en 
oeuvre en 6niettant de' fagon aussi isotrope que possible, 
20 notainment vers le fond B et depuis la source acoustique 
immergee 1, 1 ^ onde acoustique incidente s(t), de 
signature f requentielle connue. 

L' operation d ' acquisition ACQUI est, quant- a elle, mise 
25 en oeuvre en produisant, au moyen des recepteurs immerges 
21 a 24, des signaux de mesure respectifs sl(t,l), 
s2(t,l), s3(t,l) et s4(t,l) (generiquement notes Si(t)) 
en reponse a une detection d'ondes acoustiques 
respectives propagee dans I'eau, et en enregistrant ces 
30 signaux en fonction du temps t. 

Comme le montre la figure 1, la source 1 et la flute 2 
restent a distance du fond B et se trouvent 
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respectivement immergees a des profondeurs Zs et Zt , 
eventuellement egales, la flute 2 etant au moins 
approximativement alignee sur la source 1 et au moins 
approximativement parallele au fond B. 

Les profondeurs d' immersion de la source 1 et des 
recepteurs 21 a 24 sont par exemple enregistrees par des 
capteurs de pression^ la profondeur d' eau H etant par 
exemple enregistree a I'' aide d' un sondeur de bathymetrie. 

La hauteur totale de la nappe d'eau au-dessus du fond B 
etant notee H, la source 1 et la flute 2 se trouvent 
ainsi globalement, par rapport au fond B, a une hauteur 
moyenne Hm typiquement comprise entre 10 metres et 20 
metres . 

Les recepteurs 21 a 24 sont espaces les uns des autres^, 
chaque recepteur ' etant separ6 d^un recepteur voisin par 
une distance d qui peut etre constante ou Evolutive le 
20- long de la flilte 2. 

La flute 2 est consideree comme delimitee entre les deux 
recepteurs extremes 21 et 24, qui sont done equidistants 
du centre 20 de cette flute. On note 'Di.;. la distance 
25 horizontale separant l^'emetteur du i^"^® hydrophone et Zri 
la profondeur du i^^^ hydrophone. 

La phase d' analyse ANA, qui est destin^e a permettre 
1 ' identification des param^tres physiques tels que hj, 
30 pj, Cj, ofj de chaque strate Fj du fond B, exploite a cette 
fin la signature f requentielle de I'onde acoustique 
incidente s(t), et les signaux de mesure Si(t,l) des 
differents recepteurs 21 a 24. 



10 
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Le proced6 de 1' invention concentre d'abord son analyse 
sur les signaux- Si(t) propages entre la source 1 et les 
recepteurs 21 a 24. 

5 

Pour ce faire, I'onde incidente s(t);. aussi isotrope que 
possible, est emise par la source 1, de directivite 
connue, dans des fenetres temporelles d' emission 
successives disjointes, et le signal de mesure Si(t,l) de 
10 chaque recepteur est lui-meme analyse dans des fenetres 
temporelles de reception successives disjointeS;. au cours 
desquelles sont detectes les signaux propages. 

Dans ces conditions, chaque signal regu Si(t) est 
15 constitue de repliques attenuees et retardees du signal 
emis s(t) pour le tir note Tl, la structure du signal 
regu Si(t) etant ainsi fonction de la geom§trie du trajet 
multiple que I'onde acoustique parcourt dans I'eau 
(figure 5) . 

20 

Ce trajet multiple comprend principalement un trajet 
direct de I'onde acoustique depuis la source 1 jusqu'a 
chacun des recepteurs 21 a 24, un trajet dit "reflechi 
fond" de I'onde' acoustique depuis la source 1 jusqu'a 

25 chacun des recepteurs 21 a 2 4 via reflexion sur le fond, 
I'onde acoustique ayant parcouru ce trajet donnant 
naissance au signal de reflexion rfi(t,l), et un trajet 
dit "reflechi surface" de I'onde acoustique depuis la 
.source 1 jusqu'a chacun des recepteurs 21 a 24 apres 

30 reflexion a la surface de I'eau. 

Dans le cas particulier des fonds a faible hauteur d'eau 
H ou "pet its fonds", on utilisera avantageusement le 
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trajet de reflechi surface voire le trajet dit de 
"reflechi fond-surface" impliquant une reflexion de 
I'onde acoustique a la fois sur le fond et sur la surface 
entre son emission par la source 1 et sa detection par 
5 I'un des recepteurs 21 a 24. 

La mesure des temps d' arrivee respectifs de ces trajets 
multiples dans 1 ' eau permet d'une part de determiner avec 
precision la geometrie du dispositif a 1' instant de la 

10 mesure et d' autre part d' etablir des fenetres de 
reception a I'interieur desquelles se concentrera 
1' analyse des reflechis fond- 
Plus precisement, le precede de 1" invention concentre son 

15 analyse sur les signaux de mesure issus d'ondes 
reflechies une seule fois par le meme fond B. Comme 
indique sur la figure 3, doiit les differents schemas 3A a 
3C correspondent a des tirs differents^ un ordonnancement 
de plusieurs tirs cohs6cutifs permet alors de recueillir 

20 les signaux notes Rfi(t^B) issus de la reflexion de 

v.- 

1 ' onde acoustique incidente sur une meme zone de fond B 
et sous autant d' angles de rasance differents qu'il 
existe de recepteurs tels que 21 a 24. 

25 Comme le montre la figure 4, les differents signaux de 
reflexion rfi(t;.l) recueillis, lors d'un tir 1, sur 
chacun des recepteurs 21 a 24 dans une fenetre 
d' observation correspondant a la reception d'un echo 
unique sur le fond doivent done §tre respectivement 

30 associes a differents signaux de reflexion rfi(t,lMf 
rfi(t,l''), rfi(t,l'^M etc.;. recueillis, lors d'autres 
tirs 1', 1''/ l''^/' etc. sur chacun des recepteurs 21 a 
24 pour former les signaux dits de "reflechi fond" ou de 
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reflexion recales Rfi(t,B);. cette operation etant 
effectuee a une etape de recalage RECAL qui sera 
detaillee ulterieurement . 

5 En realite, chaque signal de reflechi fond Rfi(t,B) est 
lui-meme composite en ce sens qu'il est forme (figure 6) 
par la propagation, suivant un trajet multiple, de 1 ' onde 
acoustique incidente dans le fond B, donnant lieu a des 
reflexions de cette onde incidente sur les differentes 
10 strates Fj du sous-sol sub-aquatique • 

Dans ces conditions, chaque signal de reflechi fond 
Rfi(t,B) peut etre decrit sous la forme d'un ensemble de 
rayons propres notes RfSi^(t,.B)^ I'indice k se rapportant 

15 au ]c^^^^ rayon propre associe a I'^onde incidente Rfi(t,B)- 
Le rapport d' amplitude entre le signal temporel incident 
et reflechi a 1' interface Fj est directement 
proportionnel a la transf ormee de " Fourier du coefficient 
de reflexion Rj (f /Gi^( j ) ) , l''indice k se rapportant ici a 

20 I'' angle Oi^(j) avec lequel le k^^^® rayon propre vient se 
reflechir sur 1' interface F j . 

Or^ ce coefficient de reflexion Rj (f ^ 0i^ ( j ) ) dans le 
cadre d''une modelisatioh de Rayleigh^. est relie aux 
25 parametres geotechniques des differentes strates Fj~l Fj 
du sous-sol. 

Les signaux ainsi produits et exploites portent la 
signature de 1' environnement rencontre par les ondes 
30 detectees, et contiennent done une information 
suffisamment riche et pertinente du fond B pour permettre 
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d'en extraire sans ambiguite les parametres geotechniques 
des differentes strates Fj du sous-sol. 

L'echo d^une onde incidente parvenant au fond B sous un 
5 angle de rasance 9 tres rasant n'est pratiquement 
sensible qu'a I'' attenuation acoustique de la strate Fl . 

L'echo d'une onde incidente parvenant au fond B sous une 
incidence perpendiculaire au fond B n'est pratiquement 
10 sensible qu'aux changements d' impedance acoustique et 
presente notamment le defaut de ne pas pouvoir permettre 
sans ambiguite la determination du couple epaisseur- 
vitesse de chacune des strates. 

15 De fagon generale^. il est done preferable de recuperer 
des signaux ayant subi une reflexion a des angles de 
rasance assez faibles sur chacune des strates Fj 
composant le fond. 

20 Ainsi, des angles de rasance moderes sont ^necessaires 
pour decrire la premiere strate Fl et des angles de 
rasance plus importants sont necessaires pour decrire les 
couches plus profondes, F2 ^ etc, II est done apparu tres 
opportun de faire en sorte, pour 1 ' invention^- que chaque 

25 onde reflechie atteignant I'un des recepteurs 21 a 24- 
dans une fenetre temporelle de reception soit issue de la 
reflexion d'une onde incidente ayant atteint le fond B 
sous un angle de rasance 9k au plus 6gal a 70 degres et 
au moins egal a 10 degres. 

30 

Bien que les figures 5 et 6 soient a la fois totalement 
schematiques et incompletes compte tenu de leur mode de 
presentation simplifie et purement spatial, elles 
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illustrent, de fagon intuitive^ ce que sont les trajets 
multiples de I'onde incidente, respectivement dans 1 ' eau 
et dans le fond. 

5 L' importance ci-dessus discutee des angles d' emission vis 
a vis du reflechi fond peut etre observee sur les courbes 
des figures 7A a IJ, .correspondant respectivement aux 
ondes recueillies sur des recepteurs disposes a des 
distances correspondant a des angles de rasance au fond 
10 decroissant de 72 a 25 degres. 

Ainsi/ \il apparalt que la courbe qui correspond a des 

angles . de rasance forts ^ apporte peu d' information 
relativement a la courbe 7A puisque pour ces angles^ le 

15 coefficient de reflexion de chaque strate est 
pratiquement independant de 1' angle de rasance. De meme^ 
ii apparalt que les courbes 71 ■ et 7 J qui correspondent a 
des angles . de rasance f aibles n' apportent plus 
d' information puisque les ondes ne penetrent pas dans le 

20 fond. Inversement, les- figures 7B a 71 qui correspondent 
a un balayage angulaire pertinent^ present ent une forte 
variabilite permettant de retrouver sans ambigulte les 
caracteristiques du fond. 

25 De meme, alors que la signature d'une onde incidente tres 
basse frequence est difficile a maltriser^ et qu'une onde 
a haute frequence presente une capacite de penetration 
nulle dans les strates prof ondes,. il est apparu tres 
opportune pour 1* invention, d'utiliser une onde incidente 

30 presentant une composante de frequence relativement basse 
au moins egale a 100 Hz, et une composante de frequence 
relativement haute au plus egale a 8 kHz, cette onde 
incidente presentant par exemple, dans chaque fenetre 
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d' emission, une frequence continument variable entre 100 
Hz et 8 kHz. 

La relation liant la profondeur de penetration d'une onde 
5 acoustique incidente a la frequence de cette onde est 
clairement illustr6e par les figures 8A a 8J, qui 
representent , pour 5 angles de rasance 9i differents de 
1 ' onde incidente compris entre 25 et 70 degres, la 
reponse Rfi(t,B) du signal teitiporel reflechi au fond, 

10 pour un signal d' emission compris entre 400 et 3200 Hz 
(figures 8A a 8E) et compris entre 5400 et 8200 Hz 
(figure 8F a 8J) . Bien que presentant la meme resolution 
temporelle que les signaux basse frequence, les signaux 
haute frequence sont sensiblement appauvris par 

15 1' attenuation des ondes dans les di verses strates. 

Ainsi, grace a un choix approprie de 1' angle de rasance 
et grace a une variation de la frequence acoustique 
d' excitation, les parametres physiques hj, pj, Cj, cXj de 

20 chaque strate Fj peuvent §tre deduits de 1 ' onde incidente 
s (t) et des signaux de mesure Si(t,l) par 1 ' intermediaire 
d'un modele physique faisant intervenir les coefficients 
de reflexion acoustique Rj(9i^(j)) des strates Fj du fond 
B vu par les rayons propres RfSi^(t,B), leur attenuation 

25 aj, ainsi que leur dependance f requentielle . 

Autrement dit, le modele physique employe permet de tirer 
parti des dependances angulaires du . coefficient de 
reflexion et de la • dependance f requentielle de 
30 1' attenuation dans . les strates. 

En pratique, il est souvent utile de faire en sorte que 
chaque onde reflechie atteignant le centre 2 0 de la flute 
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2 dans une fenetre temporelle de reception soit issue de 
la reflexion d'une onde incidente ayant atteint le fond 
sous un angle de rasance 9i au plus egal a 70 degres et 
au moins egal a 10 degres^ 1' angle de rasance 9i pouvant 
5 meme le plus souvent etre idealement choisi dans une 
plage angulaire comprise entre 25 et 65 degres. 

Chaque recepteur 21 a 24 est alors separe de la source 
acoustique 1 par une distance D comprise entre environ 1 
10 fois et environ 4 fois la hauteur moyenne Hm, la valeur D 
= Hm concernant essentiellement le recepteur 21 le plus 
proche de la source 1^ et la valeur D = 4.Hm concernant 
plutot le recepteur 24 le plus eloigne de la source !• 

15 Dans ces conditions,, chaque recepteur, tel que 21 ^ 23, 
est espace d'un recepteur voisin, tel que 22 a 24, par 
une distance d permettant une discretisation la plus 
reguliere possible de 1' angle de rasance. 

20 En pratique, les fenetres temporelles d' emission ont 
chacune une duree au moins egale 0.1 seconde, et de 
preference plutot au moins egale a 0.5 secondes. De fagon 
generale, la duree choisie est un compromis entre un 
signal court minimisant les effets du mouvement du navire 

25 et un signal plus long augmentant le rapport signal a 
bruit. 

Ces fenetres d' emission sont repetees periodiquement avec 
une periode dependant de leur dur6e individuelle, au 
30 moins egale a 0.5 secondes, et de preference plutot au 
moins egale a 1 seconde. De fagon generale, la duree du 
silence entre deux tirs est un compromis entre une duree 
elevee permettant une s<^paration satisf aisante des 
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signaux regus et une duree faible permettant une plus 
grande resolution au sol. 

Dans le cas prefere, les signaux de mesure sont analyses 
5 par tranche de 2 secondes et traites par tir, un tir 
correspondant a un jeu de signaux forme d'un signal emis 
et des signaux de mesure correspondants regus par les 
recepteurs respectifs 21 a 24. 

10 Entre la phase d' investigation INVESTIG telle que 
precedemment decrite, et la phase d' analyse ANA qui sera 
davantage detaillee ulterieurement ^ le precede de 
1' invention comprend une phase de pre-traitement not6e 
PRE_TRAIT. 

15 

Cette phase de pre-traitement . PRE_TRAIT comprend une 
operation de filtrage FIL_COR^ une operation de calcul 
geometrique INVER_GEOM^ et une etape de recalage REGAL. 
Au . cours de 1' operation FIL_COR, chaque signal de mesure 
20 Si(t) est filtre pour normaliser ce signal,, en tenant 
compte de la signature f requentielle de 1 ' onde acoustique 
incidente s (t) . 

Au cours de 1' operation INVER_GEOM^ ia • geometrie de 
25 1' instrumentation lors du tir est recalculee et les 
fenetres temporelles de reception sont identifiees pour 
pouvoir extraire les signaux de reflexion rfi(t^l). 

Enfin, au cours de 1' etape RECAL'^ chaque signal de mesure 
30 extrait est associe a une position geographique du fond 
reperee par ses coordonnees X et et recale de-'fagon a 
tenir compte du deplacement du navire entre chaque tir et 
le tir suivant, 1" etape REGAL produisant ainsi un jeu de 
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signaux de reflexion recales Rfi(t,B) ayant vu une meitie 
zone de fond B lors de tirs differents et sous des angles 
de rasance differents. 

5 Le filtrage effectue lors de 1' operation FIL_COR, dit 
"adapte", est obtenu par estimation de l^enveloppe de 
1' intercorrelation entre le signal emis represents par 
1 ' onde acoustique incidente s(t), et chaque signal de 
mesure si (t ^ 1) . 

operation de filtrage consiste d' abord a effectuer la 
correlation : 

puis a effectuer une extraction d' enveloppe du signal 
Ci(t^l) qui se realise par exemple a I'^aide d'une phase 
de demodulation selon c^(r,/).e"'^'^'^''' , ou f c est la frequence 
centrale du signal eidis et suivie d' un filtrage passe-bas 
a la frequence de coupure Fcoup > (Fmax-Fmin) /2 pour 
obtenir le signal ei(t,l). 

Le calcul de la geometrie instrumentale exacte effectue 
lors de 1' operation INVER_GEOM 'est base sur 
1 ' identification des ondes reflechies arrivant 
successiveiuent sur les differents recepteurs 21 a 24. 

Cette identification est obtenue par exemple a partir 
d'une detection de maximum a l''interieur de fenetres 
temporelles qui sont definies a partir de la geometrie a 
priori de 1' instrumentation. 
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On compare les temps mesures (Tei^)m correspondant a cette 
detection aux temps (xei^) s issus d' un modele de rayons 
propres en milieu stratifie fournissant 1' expression de 
la reponse impulsionnelle theorique : 

5 

N 

Rie^{t)^^A,^,S{t'-Te^) Equation 1 

ou N represente le nombre de trajets acoustiques pouvant 
relier un recepteur a I'emetteur, ou t represente la 
variable temps ou xei^ represente le retard que 

10 presenter par rapport a 1' emission de 1 ' onde incidente,. 
la reception d'une onde ayant suivi le k^^^® trajet^ ou 5 
represente la fonction de DiraC;. et ou Ai^ represente 
1' amplitude relative, rapportee a celle de I'onde 
incidente emise, de I'onde ayant suivi le k^^"^® trajet et 

15 parvenant au recepteur avec un retard xei^, la notation xe 
rappelant qu' il s'agit des retards associes au trajet 
multiple de I'onde acoustique dans I'eau. 

Les retards associes au trajet direct et au r§flechi fond 
20 sont successivement donnes par : 

re. = et re, = 

^moy ^moy 

oh Cmoy est la Vitesse moyenne du son dans la colonne 
25 d' eau traversee par I'onde ayant suivi ce trajet et dans 
le cas particulier des grands fonds oti les reflexions en 
surface ne sont pas exploitees. 



Les instants theoriques prevus pour les differents rayons 
30 propres calcules sont compares aux instants d'arrivees 
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des ondes reellement detectees par les differents 
recepteurs 21 a 24. 

En cas de difference non negligeable^ une recherche de la 
5 geometrie exacte peut etre effectuee par exemple par un 
algorithme de recuit simule a 2.N variables (Di et Zri) 
en minimisant une fonction de cout exprimee a partir des 
differences entre les instants theoriques associes aux 
rayons propres et les instants de detection des ondes 
10 recueillies sur les N differents recepteurs 21 a 24. 

Dans le cas particulier des petits fonds^r d' autres 
arrivees telles que le ref lechi-surf ace et meme les 
reflechis surface-fond ou fond-surface sont detectes et 
15 permettent de s'affranchir des mesures prealables de H et 
Zs, grandeurs qui seront direct ement recherchees par 
1' algorithme d' inversion geom6trique decrit ci-dessus. 

A partir de 1 ' information alors disponible,. qui est 
20 constitute, pour chaque tir et chaque receptei^r "i", par 
le signal Si(t) issu de ce recepteur et par la geometrie 
du systeme, il convient de selectionner la plage 
temporelle correspondant a la fenetre de reception dans 
laquelle le signal regu sera analyse, c'est-a-dire la 
25 plage temporelle de reception de la seule contribution du 
signal de reflexion rfi(t,l). 

Cette selection est par exemple effectuee, dans le cas 
particulier des petits fonds, par un programme qui, pour 
30 chaque enveloppe du signal ei(t,l) issu du i^™® recepteur 
selectionne la fenetre temporelle definie par : 
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ou est un separateur symbolique place entre la borne 

temporelle inferieure et la borne temporelle superieure 
5 de la fenetre de reception selectionnee, ou Te±^ 
represente le retard associe au ref lechi . f ond^ ou xei"* 
represente le retard associe au premier reflechi surface- 
fond, ou N est un parametre dependant de la nature du 
signal emis, et ou te est la periode d' echantillonn^ge de 
10 chaque recepteur, les indices 3 et 4 etant lies a un cas 
particulier des pet its fonds dans lequel I'^ordre 
d' arrivee des rayons propres serait successivement le 
trajet , direct, le ref lechi-surf ace, le ref lechi-f ond et 
le reflechi surface-fond. 

15 

L'' operation de recalage REGAL consiste a utiliser les 
donnees de position (X,Y) avec les positions des 
instruments associe au tir correspondant de fagon a 
calculer les positions Bi"^(X,Y) du fond vu par 
20 1'' hydrophone i lors du tir 1- En repetant 1'' operation 
pour plusieurs tirs consecutifs, il devient possible de 
reconstituer le jeu de donnees Rfi(t,B) constitue des i 
reflechis fonds ayant vu le fond B(X,Y) sous les 
differents angles de rasance Gi. 

25 

A titre d'exemple, pour une configuration definie par une 
distance de 10m entre la source 1 et le premier 
hydrophone 21 de la flut^, pour une flute a 4 hydrophones 
separes les uns des autres de 10m, horizontale et situee 
30 a 10m au dessus du fond, pour un navire se deplagant a la 
Vitesse reguliere de 5m/s et avec une duree de 2 secondes 
entre tirs successifs, on appliquera successivement un 
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decalage de 1 tir pour le 2 ® hydrophone 22, de 2 tirs 
pour le 3^^""^ hydrophone 23, et de 3 tirs pour le 4^^"^^ 
hydrophone 24 coinme rappele dans le Tableau 1 suivant - 



N^'^' de tir 


Signaux 
extraits 


Decal . 


Recalage 


Tir 1 


rfi(t,l) 


0 


rfi(t,l) -> Rfi(t,B) 


Tir 1+1 


rfi(t,l+l) 


1 


rf2(t,l+l)-> Rf2(t,B) 


Tir 1+2 


rfi(t,l+2) 


2 


rf3(t,l+2)-> Rf3(t,B) 


Tir 1+3 


rfi{t,l+3) 


3 


rf4(t,l+3)-> Rf4{t,B) 



5 

Tableau 1: Principe de Top^ratioii de recalage des tirs 

Une fois identifi6e et extraite, pour chaque recepteur, 
1 ' onde detectee ayant subi une reflexion unique et par le 
10 fond B, la phase d' analyse ANA est mise en oeuvre pour 
permettre de deduire de cette onde les param^tres 
geotechniques du fond au moyen d'un modele fond6 sur 
1 " expression du coefficient de reflexion de Rayleigh. 

15 Coirane evoque precedeiriment , la ,modelisation des 
interaGtions des ondes acoustiques avec le fond et ses 
differentes strates est basee sur la representation de 
chaque strata sur une base de plusieurs parametres 
incluant I'epaisseur h de cette strate^ la densite p (par 

20 exemple en g/cm^) ^ la vitesse du son Cp/s (par exenaple en 
m/s) et 1' attenuation du son otp/s (par exemple en dB/X) ; 
les indices p et s se rapportant respectivement aux ondes 
de compression et aux ondes de cisaillement . 

25 Dans les couches sedimentaires^ la rugosite ou en'core un 
gradient de celerite peuvent dans certains cas etre 
ajoutes a ces principaux parametres physiques- 
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Partant du fait que le reflechi fond se decompose lui- 
meme coinme une sonime de rayons propres reflechis par les 
differentes strates du fond, on peut egalement en ecrire 
5 sa reponse impulsionnelle selon : 

RIb,(t) = f^4.Sit-Tb^) 



Tbi^ se rapportant au fait qu'il s'agit du trajet multiple 
10 de I'onde acoustique dans le fond. 



Considerant maintenant E(t), 1' autocorrelation de 
I'enveloppe du signal ^miS;. on obtient Rfi(t) par : 



ik=l jfc«l k=l 

I'indice s precisarit qu'il s'agit du reflechi fond 
simule . 



20 Concernant la determination de ces rayons propres, il 
s'agit en particulier de calculer les k angles de tirs 

(angle de rasance au niveau de la source* 1) des k rayons 
propres recherches. 

25 On part de 1' angle Q±^ defini par : 

e, =^tan( ^) 

On effectue alors une recherche dichotomique de-' I'' angle 
de tir Oi^(O) entre 0i° et 0i°+Gmax utilisant les 

30 relations : 
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5 



25 



cos^. (7) cos <9, (7 + 1) 

= entre les strates Fj et Fi+1, en 

notant 0i^(j) 1' angle dans la couche Fj du rayon partant 
avec un angle de tir 0i^(O)^ et 



^ ^ i2.H-Zr,^Zs) ^ y 2:NO\k).hU) ^ ^ 
' tan^/(0) tanO.'U) 

en considerant ici le rayon propre du reflechi fond 

associe k la geometrie du capteur i coitime le rayon propre 
10 se ref lechissant au plus profond sur la j^^™® strate, et 
en notant N(j,k) le nombre de reflexions du rayon k dans 
la couche Fj . 

Un test doit prealablement etre effectue pour s' assurer 
15 de 1'' existence de ce rayon N(l;.k). 

Da:ns le cas ou Ai est une fonction de la frequence, on 
a : V 

Rfs t{t,B) = f^4(t) .E(t-Tbt) 

20 k^l : 

4(t) = FFT-\4(f)) ■ ■ 

Dans le cas particulier d' interfaces planes et sans 
reflexion multiple dans aucune couche. 



wo 2005/101058 



PCT/FR2005/000587 



ou dj est la distance parcourue dans la strate Fj , 
tj=dj/Cj le temps correspondant aj 1' attenuation du son 
dans la couche exprimee en Neper/m^ et Rj(ff0i^(j)) le 
coefficient de reflexion a la j^®^® interface. 

Le coefficient de reflexion de Rayleigh associe a une 
interface f luide-f luide est defini par : 



r ( f fi\( i\\ - sm(^/(y-H)) sin(.^/(y)) 

0^-/ yj))" ri c O C 

sinC^/a + l)) sm(^/a)) 
10 ou 1 ' on utilise les celerites complexes : 



=— ^ 



1 + - • 

2.;?r.8,686 



Ce modele etant pose- en amont de toute investigation, la 
15 phase d' analyse ANA peut I'utiliser pour extraire les 
parametres geotechniques des differentes strates du sous- 
sol sub-aquatique a partir des signaux de reflexion 
recales et corriges Rfi(t,B) obtenus la fin de la phase 
PRE_TRAIT- 

20 

Pour ce faire, la phase d' analyse ANA utilise par exemple 
une technique d ' appariement de formes connue de 1 ' homme 
de metier, auquel cas elle consiste a parcourir, de fagon 
iterative, une boucle incluant une operation de 
25 modelisation MODEL, une operation de soustraction SOUST, 
et une operation d' optimisation et de retroaction 
OPTI RETR. 
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On rappelle que^ par convention, l*indice "n" sera 
utilise ci-apres pour designer de fagon generique le 
nombre courant de fois que la boucle d' iteration a ete 
parcourue depuis le debut de la phase ANA, et pour 
5 reperer toute grandeur obtenue au n^^^^ pas du parcours 
iteratif de cette boucle. 

La modelisation MODEL consiste a calculer, a chaque pas 
n, les signaux de reflexion recales fictifs respectifs 
10 Rfi^(t, Bn) que produiraient les recepteurs 21 a 24 si le 
fond B etait constitue par un fond note Bn forme de 
strates Fj definies par des paraitietres physiques fictifs 
hnjf Pnj/ Cnj, o^nj virtuelleiuent attribues a ces strates. 

15 Les paraitietres physiques fictifs initiaux hoj, Pojr Coj/ 
aoj sont par exemple choisis sur la base d' indications 
geologiques generales sur le site a explorer, ou a defaut 
au hasard dans des fourchettes de valeurs plausibles ou 
simplement possibles, les valeurs prises par les 

20 paraitietres hnj, pnj^ Cnj, cxnj pour les valeurs ^.non nulles 
de n etant definies de proche en proche par le parcours 
repete de la boucle d' iteration. 

A ces paraitietres peut etre ajoute le nombre Jn defini 
25 coitvme le nombre total de strates supposees, lors de 
1' execution du n^^^® pas de la boucle d ' iteration, avoir 
ete detectees, ce nombre evoluant jusqu'a I'obtention, a 
1' issue de 1' iteration, du nombre total exact J© de 
strates reellement d6tect6es. 

30 

Pour l^obtention des parametres exacts, la' phase 
d" analyse ANA comprend tout d^abord une operation de 
soustraction SOUST au cours de laquelle sont calcules les 
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10 



15 



20 



25 



ecarts Ani que presentent, par rapport aux signaux de 
reflexion recales fictifs RfSi^(t^ Bn) prealableitient 
calcules, les signaux de reflexion recales filtres et 



Puis, la phase d' analyse ANA comprend une operation 
d' optimisation et de retroaction OPTI_RETR a 1' issue de 
laquelle les parametres physiques fictifs sont revises en 
fonction des ecarts constates. 

Pour ce faire, le proced6 de 1' invention calcule, a 
chaque pas d' iteration n, une fonction de cout donnee par 
exemple par: 



avec Ani*= RfSi''(t, Bn) - Rfi(t,Be)- 

Cette fonction de cout est par exemple minimiaee selon un 
algorithme de recuit simule pour produire une reponse 
impulsionnelle affinee Rf Si^'^-'* (t , Bn+i) definissant un 
nouveau fond Bn+i utilis6 dans un nouveau pas d' iteration 
de la phase ANA, le processus convergeant vers la 
production de parametres exacts hej, Pejf Cej, cXej des 
dif f erentes strates j . 

Le tableau ci-dessous donne un exemple de cas concret 
traite - 



corriges Rfi(t,Be) . 
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Geometrie : 






H=25m^ Zs=7m, 


Zri=8m, Di=12m, di 


=5in, Je = 5 


Fond : 






Cpl^=1750ni/s; 


pl^=1.8; al°=0.5; 


hl''=0.75iti; 


Pr^9®— 1 Q7nm / o • 
L^pZ — ±^ / \JUX/ 3 f 


pz — x.yj; cxz — 


Liz —X . / DUXf 


Cp3^=17 93m/s; 


p3^=1.85; a3^=0.5; 


h3^=0. 4 5ni; 


Cp4^=2050m/s; 


p4^=2.; a4^=0.38 


; h4^=2.25m; 


Cp5^=3000m/s; 


p5^=2.3; a5^=0.2; 


Cs5^=1300m/s; a5s^=0.3 



Coirime le Gomprendra I'homme du metier, au lieu d'utiliser 
une fonction de cout telle que K(n) calculee au moyen une 
distance euclidienne^ la phase d' analyse ANA pourrait 
5 tout aussi bien utiliser la recherche d'un maximum de 
correlation entre les signaux Rfsi^(t;r Bn) et Rfi(t,Be). 

Comme neanmoins il s'agit^ dans I'un et 1' autre de ces 
deux cas, de reveler des ecarts entre ces signaux et de 

10 faire evoluer le signal RfSi^(t, Bn) dans un sens propre a 
reduire ces ecarts, les operations pr6cedemment evoquees 
de soustraction SOUST et de formation de tels ecarts 
doivent etre ici comprises dans un sens fonctionnel large 
incluant notaitiment une telle recherche de maximum de 

15 correlation. 
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REVENDICATIONS. 

1. Precede de caracterisation geotechnique d'un fond 
5 sub-aquatique (B) , tel qu'un fond sous-marin, recouvert 
d'une nappe d'eau de hauteur totale (H) determinee, et 
comprenant une pluralite de strates (Fl, F2) formant 
entre elles des interfaces de separation^ presentant 
chacune des parametres physiques propres (hj, pj^ Cj, aj) 
10 et s'etendant a des profondeura di verses en dessous d'une 
premiere strate (Fl) qui forme une interface de 
separation avec l^eau, ce precede comprenant au moins une 
phase d' investigation (INVESTIG), une phase de pre- 
traitement (PRE_TRAIT) et une phase d'analyse (ANA), la 
15 phase d' investigation (INVESTIG) incluant elle-m§me une 
operation (EXCIT) d' excitation acoustique;. mise en oeuvre 
en emettant au moins vers le fond (B) depuis une source 
acoustique (1) immergee, une onde acoustique incidente 
(s (t) ) de signature frequentielle connue, une operation 
20 d' acquisition (ACQUI) mise en oeuvre en produisant, au 
moyen d'une flute (2) d'au moins quatre recepteurs >{21 a 
24) immerges, des signaux de mesure respectifs (Si(t,l)) 
resultant d'une detection d'ondes acoustiques respectives 
reflechies par le fond (B) , et une operation de releve 
25 (RELEV) de donnees geometriques (Zs (t) , Zr(t), H(t), 
X(t), Y(t)), la flute (2) etant au . moins 
approximativement alignee sur la source (1) et dista-nte 
du fond (B), les recepteurs (21 a 24) etant espaces les 
uns des autres, et I'onde incidente (s(t)) etant emise 
30 dans chaque fenetre d'une succession de fenetres 
temporelles d' emission disjointes, et presentant une 
variation temporelle de frequence a I'interieur de chaque 
fenetre d' emission, caracterise en ce que la phase 
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d' investigation (INVESTIG) inclut une operation de 
deplacement (DEPL) , concomitante avec les operations 
d'excitation (EXCIT) et d ' acquisition (ACQUI) , et mise en 
oeuvre en deplagant simultanement la source acoustique (1) 
5 et la flute (2) dans la nappe d'eau^ en ce que la phase 
de pr6-traitement (PRE_TRAIT) est mise en ceuvre en 
deduisant, a partir des signaux de mesure Si(t,l), de la 
signature de 1 ' onde acoustique incidente {s(t)), et des 
donnees geometriques (Zs (t) , Zr(t), H{t), X(t), Y(t)), 
10 des signaux de reflexion rfi(t,l) ayant atteint des 
recepteurs differents (21 a 24) a la suite d'une m^me 
emission d'onde incidente, et en regroupant, sous forme 
de series de signaux de reflexion recales (Rfi(t,B)), des 
signaux de reflexion rfi(t,l) ayant atteint, sous des 
15 angles de rasance {Q±) differents, des recepteurs 
differents (21 a 24) en provenance d'une meme zone du 
fond a la suite de differentes emissions de I'onde 
incidente, et en ce que la phase d' analyse (ANA) est mise 
,en CBuvre en d^duisant les parametres physiques (hj, pj, 
20 Cj, aj) de chaque strate (Fl a F4) du fond, a,.partir des 
series de signaux de reflexion recales (Rfi(t,B)), par 
exploitation de la dependance angulaire et f requentielle 
' du coefficient de reflexion de I'onde incidente sur 
chaque interface, chaque signal de reflexion recale 
25 (Rfi(t,B)) etant interprets comme une somme d'arrivees 
retardees et attenuees de rayons propres correspondant a 
la reflexion de I'onde incidente par les differentes 
strates de la zone du fond a laquelle correspond ce 
signal de reflexion recale (Rfi(t,B)). 

30 

2. Procede de caracterisation geotechnique suivant la 
revendication 1, caracterise en ce que chaque onde 
ref lechie atteignant un recepteur est issue de la 
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10 



15 



20 



25 



reflexion d'une onde incidente ayaht atteint le fond sous 
un angle de rasance (9i) au plus egal a 70 degres et de 
preference au plus egal a 65 degres. 

3. Procede de caracterisation geotechnique suivant 
I'une quelconque des revendications precedentes, 
caracterise en ce que chaque onde reflechie atteignant un 
recepteur est issue de la reflexion d'une onde incidente 
ayant atteint le fond sous un angle de rasance {e±) au 
moins egal a 10 degres et de preference au moins egal a 
25 degres. 

4. Procede de caracterisation geotechnique suivant 
I'une quelconque des revendications precedentes, 
caracterise en ce que chaque recepteur (21 a 24) est 
separe d'un recepteur voisin (22 k 24) par une distance 
(d) correspondant a une difference d'' angle de rasance au 

moins 6gale ^ 2 degres, et de preference au moins egale a 
5 degres, des signaux de reflexion respectivement reQus 
par ce recepteur et par le recepteur voisin. 

5. Procede de caracterisation geotechnique suivant 
I'une quelconque des revendications precedentes, 
caracterise en ce que 1 ' onde incidente presente, dans 
chaque fenetre d' emission, au moins une composante de 
frequence relativement basse, dont la frequence est au 
moins egale a 100 Hz. 



6. Procede de caracterisation geotechnique suivant 
30 I'une quelconque des revendications pr6cedentes, 
caracterise en ce que I'onde incidente pr6sente, dans 
chaque fenetre d' emission, au moins une composante de 
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frequence relativement haute, dont la frequence est au 
plus egale a 8 kHz. 

7. Precede de caracterisation geotechnique suivant 
5 I'une quelconque des revendications precedentes^ 
caracterise en ce que 1 ' onde incidente presente, dans 
chaque fenetre d' emission, une frequence continument 
variable entre une composante de frequence relativement 
basse et une composante de frequence relativement haute. 

10 

8- Precede de caracterisation geotechnique suivant 
I'une quelconque des revendications . precedentes, 
caracterise en ce que les fenetres temporelles d' emission 
ont chacune une duree au moins egale a 0-1 seconde et de 
15 preference au moins egale a 0.5 seconde . 

9. Procede de caracterisation geotechnique suivant 
I'une quelconque des revendications pr6cedentes, 
caracterise en ce que les fenetres temporelles d' emission 

20 sont repetees periodiquement avec une perio4e au moins 
egale a 0.5. seconde et de preference au moins egale a 1 
seconde - 

10. Procede de caracterisation geotechnique suivant 
25 I'une quelconque des revendications precedentes, 

caracterise en ce que les parametres physiques sont au 
moins en partie choisis dans 1' ensemble de parametres 
comprenant I'epaisseur (h) , la densite (p) , la vitesse du 
son (Cp, Cs) , I'attenuation du son (ap, as), la rugosite, 
30 et un gradient de celerite du son. 

11. Precede de caracterisation geotechnique suivant 
I'une quelconque des revendications precedentes. 
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caracterise en ce qu'il comprend, en araont de la phase 
d' analyse (ANA), une operation (CAPT) de determination de 
la profondeur d' immersion (Zs) de la source acoustique 
et/ou de la profondeur d' immersion (Zr) d^un point de la 
5 flute et/ou de la hauteur d''eau H. 

12. Precede de caracterisation geotechnique suivant la 
revendication 10, caracterise en ce que la phase 
d' analyse (ANA) est mise en oeuvre en parcourant, de fagon 

10 iterative, une boucle incluant une operation de 
modelisation (MODEL_FOND) au cours de laquelle sont 
calcules, a partir de parametres physiques fictifs (hnj, 
Pnj/ Cnj/ anj) virtuellement attribues a la premiere strate 
* (Fl) au moins de la plurality de strates, et pour les 

15 differents recepteurs (21 a 24), des signaux de reflexion 
recales fictifs respectifs (Rfsi''(t, Bn) ) issus de signaux 
de mesure virtuellement produits par ces recepteurs (21 a 
24) en presence de ces parametres fictifs (hnj^ pn^r Cnj, 
anj), une operation de soustraction (SOUST) au cours de 

20 laquelle sont formes des* ecarts (Ani) constate^ entre les 
signaux de reflexion recales respectifs (Rfi(t,B)) issus 
des signaux de mesure (Si(t, 1)) reellement produits par 
ces recepteurs (21 a 24) et les signaux de reflexion 
recales fictifs correspondants (Rfsi''(t, 'Bn)), et une 

25 operation d ' optimisation et de retroaction (OPTI_RETR) a 
1' issue de laquelle les parametres physiques fictifs sont 
revises en fonction des ecarts constates. 



wo 2005/101058 



PCT/FR2005/000587 



1/6 




wo 2005/101058 



PCT/FR2005/000587 



2/6 



INVESTIG 



DEPL 




EXCIT 




ACQUI 




CAPT 




REFER 



RELEV 



S(t) 



Si(t,l) 



Zs(t) 
Zr(t) 
H(t) 



X(t) 
Y(t) 



V 



FIL COR 



i 



e j(tl) 



V 



INVERT GEOM 



rfi(U) 



V 



Zs.Zri 
H,Di 



REGAL 



PRE TRAIT 



V 



Rfi(t,B) 



ANA 



n=0 



Bn+1 




OPT! RETR 



Bn 



V 



MODEL FOND 



R fi"(t,Bn 



SOUST 



mi 



V 



Fig. 2 



hej(X,Y),pej(X,Y),Cej(X,Y),aej(X,Y) 



wo 2005/101058 



PCT/FR2005/000587 



3/6 




i: \ /Si(ti.ii) 




Fig. 3A 



21 



B(X,Y) 
22 23 



B(X,Y) 



Fig.SB 



24 




Fig. 3C 



rf^ai) rfgCtl) rfgCt,!) rf4(t,l) 




Fig.4A 



Rf2(t,B) Rf3(t,B) 




Rf4(t.B) 



Fig. 4B 



wo 2005/101058 



PCT/FR2005/000587 



4/6 




A 



wo 2005/101058 



PCT/FR2005/000587 



5/6 





wo 2005/101058 



PCT/FR2005/000587 




INTERNATIONAL SEARCH REPORT 



intei 



rr^^ial Application No 

PCT/FR2005/000587 



A. CLASSIRCATION OF SUBJECT MATTER 

IPC 7 601V1/38 



According to International Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searclied (classification system followed by classification symbols) 

IPC 7 GOIV 



Documentation searclied other tlian minimum documentation to the extent that such documents are included In the fields searched 



Electronic data base consulted during the international search (name of data base and, where practical, search terms used) 

EPO-Internal , WPI Data 



0. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ° Citation of document, with indication, where appropriate, of the reievant passages 



Relevant to claim No. 



us 6 151 556 A (ALLEN KENNETH PAUL) 
21 November 2000 (2000-11-21) 
cited in the application 
column 3, line 61 - column 4, 
column 6, line 13 - column 8, 
claim 1; figures 1,3 



1-12 



line 20 
line 40 



EP 0 215 703 A (INST FRANCAIS DU PETROL) 

25 March 1987 (1987-03-25) 

cited in the application 

page 3, line 5 - page 5, line 15 



WO 01/90782 A (VINOE VETLE ;NORSAR (NO)) 
29 November 2001 (2001-11-29) 
page 6, line 5 - page 7, line 9 
claim 1 



1-12 



1-12 



-/-- 



Further documents are listed in the continuation of box C. 



Patent family members are iisted in annex. 



** Special categories of cited documents : 

'A' document defining the general state of the art which is not 
considered to be of particular relevance 

'E' earlier document but published on or after the International 

filing date 

■L' document which may throw doubts on priority claim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use, exhibition or 
other means 

■P' document published prior to the international filing date but 
later than the priority date claimed 



■T* later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
Invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an Inventive step when the document is tai<en alone 

■y* document of particular relevance; the claimed Invention 

cannot be considered to Involve an Inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

document member of the same patent family 



Date of the actual completion of the international search 



14 July 2005 



Date of mailing of the international search report 



02/08/2005 



Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan 2 
NL - 2280 HV Rijswijk 
Tel. (+31-70) 340-2040, Tx. 31 651 epo nl, 
Fax: (+31-70) 340-3016 



Authorized officer 



Schneiderbauer, K 



Form PCT/ISA/210 (second sheet) (January 2004) 



page 1 of 2 



INTERNATIONAL SEARCH REPORT 



Interi 



I No 



n^^iai Application I 

PCT/FR2005/000587 



C.(Continuation) DOCUMENTS CONSIDERED TO BE RELEVANT 



Category " Citation of document, with indication, wliere appropriate, of tiie relevant passages 



Relevant to claim No. 



us 4 357 689 A (RUEHLE WILLIAM H) 

2 November 1982 (1982-11-02) 

column 3, line 65 - column 5, line 41 



1-12 



Form PCT/ISA/210 (continuation of second sheet) (January 2004) 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 

iiiiormation on patent family members 


Intern^^al Application No 

PCT/FR2005/000587 


Patent document 
cited in searcli report 


Publication 
date 


Patent family 
member(s) 


Publication 
date 



US 


6151556 


A 


21-11-2000 


AU 


772262 


B2 


22-04-2004 










AU 


5743000 


A 


09-01-2001 










Di\ 


nm 1 77^^ 

UU 1 i / / 0 


A 
rt 
















Mi 












FA 


O / DU 


D J. 












FP 




MX 


1 OR— 9009 














A 
rt 


1 O— 09— 9009 










wu 


nn7Q'^n'^ 

uu / youo 


A1 


90_i 9—9000 


EP 


0215703 


A 


25-03-1987 


FR 


2586820 


Al 


06-03-1987 










DE 


~" -3667333 


Dl 


11-01-1990 










EP 


0215703 


Al 


25-03-1987 










JP 


5058514 


B 


26-08-1993 










JP 


62098286 


A 


07-05-1987 










NO 


863439 


A ,B, 


02-03-1987 










US 


4796238 


A 


03-01-1989 



WO 0190782 A 29-11-2001 



NO 


20002712 


A 


26- 


-02- 


-2002 


AU 


5895401 


A 


03- 


-12- 


-2001 


OA 


2410272 


Al 


29- 


-11- 


-2001 


EP 


1301807 


Al 


16- 


-04- 


-2003 


WO 


0190782 


Al 


29- 


-11- 


-2001 


US 


2002010549 


Al«,. 


24- 


-01- 


-2002 



US 4357689 A 02-11-1982 



US 


4146870 


A 


27-03-1979 


AU 


510036 


B2 


05-06-1980 


AU 


2736077 


A 


01-02-1979 


BR 


7704936 


A 


25-04-i978 


. CA 


1096963 


Al 


03-03-1981 


DE 


2734091 


Al 


02-02-1978 


DK 


338777 


A 


29-01-1978 


FR 


2360087 


Al 


24-02-1978 


GB 


1583042 


A 


21-01-1981 


JP 


53019901 


A 


23-02-1978 


NL 


7708214 


A 


31-01-1978 


NO 


772567 


A ,B, 


31-01-1978 


NO 


821289 


A 


31-01-1978 


NZ 


184749 


A 


19-10-1981 



Form PCT/ISA/210 (patent family annex) (January 2004) 



RAPPORT DE RECHERCHE INTERNATIONALE 



A. CLASSEMENT DE L'OBJET OE LA DEIUIANDE 

CIB 7 GOlVl/38 



Dema^^internationale No 

PCT/FR2005/000587 



Selon la classification Internationale des brevets (CIB) ou a lafols selon ia classifica tion nationale et la CIB 
B. DOMAINES SUR LESQUELS LA RECHERCHE A PORTE 

Documentation minimale consultee (systeme de classification suivi des symboles de classement) 

CIB 7 GOIV 

Documentation consultee autre que la documentation minimale dans la mesure oCj ces documents relevent des domalnes sur lesquels a porte la recherche 
Base de donnees electronique consultee au cours de ia recherche Internationale (nom de la base de donnees. et si realisable, termes de recherche utilises) 

EPO-Internal , WPI Data 



C. DOCUMENTS CONSIDERES COMME PERTINENTS 



Cat6gorie ° 


Identification des documents cites, avec, le cas eciieanl, I'indicatlon des passages pertinents 


no. des revendications vis^es 


A 


US 6 151 556 A (ALLEN KENNETH PAUL) 
21 novembre 2000 (2000-11-21) 
cite dans la demande 

colonne 3, ligne 61 - colonne 4, Hgne 20 
colonne 6, ligne 13 - colonne 8, ligne 40 
revendication 1; figures 1,3 


1-12 


A 


EP 0 215 703 A (INST FRANCAIS DU PETROL) 

25 mars 1987 (1987-03-25) 

cite dans la demande 

page 3, ligne 5 - page 5, ligne 15 


1-12 


A 


WO 01/90782 A (VINJE VETLE ;NORSAR (NO)) 
29 novembre 2001 (2001-11-29) 
page 6, ligne 5 - page 7, ligne 9 
revendication 1 

-/-- 


1-12 


X Voir la suite du cadre C pour la fin de la liste des documents )( Les documents de families de brevets sent Indiqu6s en annexe 



' Categories speciales de documents cites: 

A' document defrnissant l*etat g6n6ral de la technique, non 
considere comme partlculierement pertinent 

E' document anterieur, mais publl6 k\a dale de d6p6t international 
ou apres cette date 

L" document pouvant jeter un doute sur une revendication de 
priority ou cite pour determiner la date de publication d'une 
autre citation ou pour une raison speciale (telle qu'lndiquee) 

O' document se referant a une divulgation orale, a un usage, h 
une exposition ou tous autres moyens 

P" document publie avant la date de dep6t International, mais 
posterieurement a la date de priorite revendiquee 



document ulterieur public apres la date de dep6t international ou la 
date de priorite et n'apparlenenant pas ^ I'^tat de la 
technique pertinent, mais cite pour comprendre le principe 
ou la theorie constltuant la base de I'invention 

document partlculierement pertinent; I'inven tion revendiquee ne peut 
6lre consld6r6e comme nouvelle ou comme Impllquant une actlvlte 
Inventive par rapport au document consld6r6 isol6ment 

document partlculierement pertinent; I'inven tion revendiquee 
ne peut fetre consideree comme impllquant une actlvlte inventive 
iorsque le document est assocle a un ou plusleurs autres 
documents de m§me nature, cette combinalson 6tanl evidente 
pour une personne du metier 

document qui fait partie de la me me famille de brevets 



Date & laquelle la recherche Internationale a 6le effectivement achevee 

14 juillet 2005 


Date d'expedrtlon du present rapport de recherche Internationale 

02/08/2005 


Nom et adresse postale de I'adminlstration chargee de la recherche Intematlonale 
Office Europeen des Brevets, P.B. 5818 Patentlaan 2 

NL-2280 HV Rijswijk 

Tel. (+31-70) 340-2040, Tx. 31 651 epo nl, 

Fax: (+31-70) 340-3016 


Fonctionnaire autorls^ 

Schneiderbauer , K 



page 1 de 2 



RAPPORT DE RECHERCHE INTERNATIONALE 



C.(suite) DOCUMENTS CONSIDERES COMME PERTINENTS 



Dema^pnternationale No 

PCT/FR2005/000587 



Categorle Identification des documents cites, avec, le cas echeant, I'indication des passages pertinents 



no. des revendications visees 



us 4 357 689 A (RUEHLE WILLIAM H) 

2 novembre 1982 (1982-11-02) 

colonne 3, ligne 65 - colonne 5, Hgne 41 



1-12 



Formulalre PCT/ISA/210 (suite de la deuxifeme feuille) (Janvier 2004) 



page 2 de 2 



RAPPORT DE RECHERCHE INTERNATIONALE 

Renseignements relatifs aUA membres de families de brevets 



Dema^plnternationaie No 

PCT/FR2005/000587 



Document brevet cite 
au rapport de recherche 



Date de 
publication 



l\/lembre(s) de la 
famine de brevet(s) 



Date de 
publication 



us 6151556 



21-11-2000 



AU 772262 B2 

AU 5743000 A 

BR 0011773 A 

CA 2375495 Al 

EA 3760 81 

EP 1204885 Al 

NO 20016144 A 

WO 0079303 Al 



22-04- 
09-01- 
02-04- 
28-12- 
28-08- 
15-05- 
19-02- 
28-12- 



-2004 
-2001 
-2002 
-2000 
-2003 
-2002 
-2002 
-2000 



EP 0215703 



25-03-1987 



FR 
DE 
EP 
JP 
JP 
NO 
US 



2586820 Al 

3667333 Dl 
0215703 Al 
5058514 B 
62098286 
863439 
4796238 



A 
A 
A 



,B, 



06- 03-1987 
11-01-1990 

25- 03-1987 

26- 08-1993 

07- 05-1987 

02- 03-1987 

03- 01-1989 



WO 0190782 



29-11-2001 



NO 

AU 
CA 
EP 
WO 
US 



20002712 A 
5895401 A 
2410272 Al 
1301807 Al 
0190782 Al 
2002010549 Al 



26-02-2002 
03-12-2001 
29-11-2001 
16-04-2003 
29-11-2001 
24-01-2002 



US 4357689 



A 



02-11-1982 



US 
AU 
AU 
BR 
CA 
DE 
DK 
FR 
GB 
JP 
NL 
NO 
NO 
NZ 



A 
A 



4146870 A 
510036 B2 

2736077 

7704936 

1096963 Al 

2734091 Al 
338777 A 

2360087 Al 

1583042 A 
53019901 

7708214 
772567 A 
821289 A 
184749 A 



A 

A 



27-03- 
05-06- 

01- 02- 
25-04- 
03-03- 

02- 02- 
29-01- 
24-02- 
21-01- 
23-02- 
31-01- 
31-01- 
31-01- 
19-10- 



1979 
1980 
1979 
1978 
1981 
1978 
1978 
1978 
1981 
1978 
1978 
1978 
1978 
1981 



Formulaire PCT/ISA/210 (annexe families de brevets) (Janvier 2004) 



